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We present a new formulation of the theory of an arbitrary quantum impurity in an antiferro- 
magnet, using the 0(3) non-linear sigma model. We obtain the low temperature expansion for 
the impurity spin susceptibilities of antiferromagnets with magnetic long-range order in the ground 
state. We also consider the bulk quantum phase transition in d = 2 to the gapped paramagnet (d 
is the spatial dimension): the impurity is described solely by a topological Berry phase term which 
is an exactly marginal perturbation to the critical theory. The physical properties of the quantum 
impurity near criticality are obtained by an expansion in (d — 1). 



I. INTRODUCTION 

Recent papers 1 ^ have presented a general field theo- 
retical discussion of the low energy properties of a spin 
S impurity embedded in an antiferromagnet or a super- 
conductor which is in the vicinity of a bulk spin-ordering 
quantum transition. These studies were motivated by a 
variety of recent experiments studying Zn and Ni impu- 
rities in the cuprate superconductors and spin-gap com- 
pounds. The motivations and prior work have been dis- 
cussed in some detail in Ref. (hereafter referred to as I), 
and so will not be repeated here. Further theoretical 3,4 , 
numerical^S*i, and experimental^ work on these issues 
has also appeared, and we will discuss some of these re- 
sults below. There has also been related work on impu- 
rity models in systems with fermionic excitation o 9 ' 10 . 

The purpose of this paper is to provide additional re- 
sults for the same quantum impurity problem using a dif- 
ferent field-theoretic formulation. The results in I were 
obtained using an expansion in (3 — d), where d is the 
spatial dimensionality. Stimulated mainly by the recent 
results of Hoglund and Sandvik-, we have succeeded in 
obtaining a formulation which permits an expansion in 
e = d — 1, and this will be described in the present pa- 
per. The universal scaling structure we shall describe 
below in the (d — 1) expansion turns out to be identical 
to that obtained in I using the (3 — d) expansion. This is 
strong evidence that a fixed point with the same scaling 
properties does indeed describe the physical situation in 
d=2. 

Throughout this paper, we will implicitly assume in 
our discussion that 1 < d < 3, unless stated otherwise. 
The only exception is Appendix^ where we will present 
results in d — 3. Also we will set e = (d — 1), whereas in 
I the same symbol was used for (3 — d) . 

Let us outline the main results of I and those that will 
be presented here. Consider a simple two-dimensional 
quantum antiferromagnet which undergoes a quantum 
transition from a magnet Neel state to a gapped, con- 
fining paramagnet with only integer spin excitations e.g. 



a model of coupled spin ladders^. We tune the antifer- 
romagnet across this transition with a generalized cou- 
pling go, such that there is Neel order for go < g c , and 
a gapped paramagnet for go > g c . Insert an arbitrary 
quantum impurity (e.g. a vacancy) which leads to a net 
deficit or excess of spin S in its vicinity (after account- 
ing for the sublattice alternation). At a temperature T 
above the gapped paramagnet phase, with go > g c and 
a spin gap A, this impurity will contribute an impurity 
spin susceptibility 
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with exponentially small corrections as T — > (we set 
h = kg = 1 and have absorbed factors of the gyromag- 
netic ratio and the Bohr magneton in the definition of the 
external magnetic field) . We can view 1)1. lfl as a definition 
of the value of S (which must be an integer or half-odd- 
integer) for the quantum impurity. In the magnetically 
ordered phase with go < g c , there are much stronger cor- 
rections to the isolated impurity behavior because of the 
presence of broken spin rotation symmetry at T = and 
gapless excitations in the bulk; in dimensions d < 2 the 
symmetry is restored at any T > 0, and corrections to 
the impurity susceptibility can be written in the scaling 
form 2 
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where p s is the spin stiffness of the bulk ordered anti- 
ferromagnet in the absence of the impurity, and c is the 
bulk spin- wave velocity. In the limit T — > 0, it was ar- 
gued in I that $(0) = S 2 /3 exactly. This prediction has 
been verified recently in the numerical study by Hoglund 
and SandvikSi. On the basis of the (3 — d) expansion, the 
subleading behavior <&(y — ► 0) = S 2 /3 + C3IJ, with C3 a 
universal number, was proposed for d = 2 in I. Hoglund 
and Sandvik^ also tested this subleading behavior, and 
argued that it did not hold — instead they proposed the 
presence of ln(l/T) term. We will show here that their 
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proposal is indeed correct, and that precisely in d = 2, 
the behavior in the limit T <C p s is 



that 
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with C3.4 unknown universal constants. The ln(l/y) de- 
pendence is special to d = 2 and does not appear at any 
finite order in the (3 — d) expansion, and this is the rea- 
son it was overlooked in I. Subleading singularities in the 
small y expansion do appear naturally in the (d — 1) ex- 
pansion presented in this paper. We also note here that 
(|1.3|) was obtained with no assumptions on the value of 
S: the S dependencies in the co-efficients are therefore 
exact. 

The subdominant ln(l/T) dependence implied by (|1.3|) 
(and the anomalous powers of y in 113.1411 1 is a conse- 
quence of spin- wave Goldstone fluctuations in 1 < d < 2, 
and does not involve the critical singularities at go = g c 
in an essential way. Consequently, in d — 2, this ln(l/T) 
dependence should also be present in antifcrromagnets 
with go <C g c , which arc not especially close to any quan- 
tum critical point. In this situation can surmise that 
(|1.3|) implies 
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where, in general, the constants Ci,2 are non-universal; 
only as we approach the quantum critical point and 
p s — > do C1.2 become universal, and then (|1.4fl is seen to 
be consistent with l|1.3(l . The ln(l/T) correction in (|1.4|) 
is related to the logarithmic frequency dependencies dis- 
cussed by Nagaosa et at— and Chernyshev et alm&. To 
the extent that sharp spin-waves are also present in or- 
dered metallic antiferromagnets, (|1.4fl may also apply to 
such systems^. 

As we will see shortly, <|1.4H is obtained for the case 
where the coupling between the impurity and the bulk 
antiferromagnet has scaled to infinity. This implies that 
at low energies the impurity moment is effectively locked 
along the direction of the local orientation of the bulk an- 
tiferromagnetic order. While such locking is appropriate 
near the quantum critical point, it is not a priori clear 
whether it should also hold at low T above a well ordered 
antiferromagnet with go <C g c . We will briefly address 
this issue by also examining the case of finite coupling 
(see Appendix [BJ : we find that the co-efficient of the 
(1/ p s ) ln(l/T) term in (|1.4|l remains universal, but there 
are non-universal corrections to the Tln(l/T) term. 

A separate category of our results concern Ximp at the 
quantum critical point, go = g c . These correspond to the 
large y, T> p s , limit of (|1.2fl . Here, it was argued in I 
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with C\ a universal number. A (3 — d) expansion for C\ 
was provided in I, and it contained non-trivial corrections 
to the free moment value of S(S + l)/3. Sushko^ has 
questioned the existence of such corrections, but we re- 
ply to his arguments in Appendix IDl The present paper 
will show that l|1.5H is obeyed also in the (d — 1) expan- 
sion: in this case the (d— 1) expansion provides terms as 
corrections to the 'classical' moment value of S 2 /3, and 
details of this appear in the body of the text, and the 
final result is in (|4.13|) . 

A number of other results for universal properties of 
the impurity correlations were provided in I using the 
(3 — d) expansion. All of these can also be computed in 
the (d— 1) expansion, and in every case we find complete 
agreement in the structure of the scaling properties. De- 
tails of such computations also appear in the body of this 
paper. 

The following section will introduce the non-linear 
sigma model field theory which describes the dynamics of 
an impurity in a quantum antiferromagnet. Section IlIII 
will then discuss the perturbative structure of this the- 
ory, with details of the perturbative computations ap- 
pearing in Appendix [S] We will show how to deduce low 
temperature properties using this perturbation theory. 
Finally, Section IIVI presents a renormalization analysis 
which allows us to deduce the physical characteristics of 
the critical point. 



II. FIELD THEORY 



This section will introduce the field-theoretical formu- 
lation of the quantum impurity dynamics which enables 
an expansion of its universal properties in the (d — 1) 
expansion. In contrast to our earlier (3 — d) expansion, 
which used a 'soft-spin' formulation of the bulk antifer- 
romagnetic fluctuations, the present (d — 1) expansion 
will use the 'fixed- length' representation of the 0(3) non- 
linear sigma model. 

We begin by recalling our earlier 'soft-spin' formula- 
tion. The bulk spin fluctuations of the antiferromag- 
net are represented by the real field 4> a (x,T), with a = 
1 ... 3 an index representing the spin component, x a d- 
dimensional spatial co-ordinate, and r is imaginary time. 
The impurity spin is placed at the origin of co-ordinates 
x = 0, and is represented by a unit length field n Q (r), 
and the bulk and impurity fluctuations are coupled in the 
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The transition in the bulk antiferromagnet is described 
by the usual 0* theory which is represented by Sb [4> a ] 
as in I. The first term in the impurity action 5; mp is the 
Berry phase of the impurity at site r: and A a (n) is a 
'Dirac monopole' function which satisfies 
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Finally, 7 is the coupling between the impurity and bulk 
degrees of freedom which will be important in our con- 
siderations below. At the 7 = fixed point, the bulk 
and boundary degrees of freedom are decoupled, and the 
coupling 7 is a relevant perturbation with scaling dimen- 
sion (3 — d — 7])/ 2 (j) is the anomalous dimension of the 
bulk critical point, and its value is very close to zero). 
The small scaling dimension of 7 near d = 3 was the key 
feature which was used to generate the (3 — d) expansion 
of the coupled bulk-impurity theory. 

Let us now turn to spatial dimensions just above d = 1. 
For the bulk theory, it is known that an expansion of 
the critical properties can be generated in a e = d — 1 
expansion by representing the bulk spin fluctuations by 
a fixed-length field N a (x,r) cx (f) a {x, r), and with the 
action of the 0(3) non- linear sigma modeli^. At the same 
time, the coupling 7 has a scaling dimension w 1, and 
so is strongly relevant near 7 = 0. This suggests that 
a better approach now would be to start near the 7 = 
00 limit. At 7 = 00, the impurity degrees of freedom 
n Q (r) would follow the bulk spin fluctuations perfectly, 
and hence n a (j) = N a (x = 0,r). In this manner we 
obtain the central field theory of interest in this paper 

Z = J VN a (x, t)S (N* - I) exp (S h [N a ] - S imp ) 



S h [N a ] = 



,1/T 

, d d x / dr 
2cg J J 

1/T 



(d T N a y + c 2 (v x N a y 



d.T 



iSA a (n) 



dn a (r) 
dr 



with n a (r) = N a (x = 0, r) 



(2.3) 



We will set c = 1 in the remainder of the paper as it 
does not appear in any essential manner in any of our 
expressions, and it can be easily re-inserted by dimen- 
sional analysis. The Berry phase in S lmp is invariant un- 
der global spin rotations and is independent of the gauge 



choice for A a . Using an analysis very similar to that pre- 
sented in I, it can be shown, order by order in (d— 1), that 
there are no relevant perturbations to the terms shown 
in l|2.3l) at the quantum critical point. Furthermore, the 
Berry phase <Si mp turns out to be an exactly marginal per- 
turbation to the bulk critical point, whose coupling con- 
stant (S) is protected by its topological nature. There 
is only a single remaining coupling constant in Z, and 
that is the bulk coupling go, and its renormalization is 
unaffected by the presence of a single impurity spin. As 
in Ref. Il5l all bulk and impurity spin correlations can be 
computed order by order in go in a diagrammatic per- 
turbation theory. We defer discussion of the structure 
of this diagrammatic expansion to Appendix [S] We note 
here that this perturbation theory makes no assumptions 
on the value of the impurity spin S, and the Berry phase 
is fully accounted for at each order in the perturbation 
theory in g . 

It is worth noting here that a perturbation theory in 
powers of go can also be generated for an arbitrary value 
of 7, with n a (r) 7^ N a (x — 0, r) (no expansion in 7 or 
S is needed here). In this case there is an additional 
gapped excitation corresponding to the deviation of the 
impurity spin from the bulk antiferromagnetic spin fluc- 
tuations (the gap of this excitation is of order 7). This 
perturbation theory is somewhat more cumbersome and 
is discussed briefly in Appendix 151 



III. PERTURBATION THEORY AT LOW T 

Before embarking upon the subtleties of a renormal- 
ization group analysis (and the associated analytic con- 
tinuation in dimensionality), it is useful to examine the 
expressions in Appendix IA II directly in 1 < d < 2, in a 
regime where perturbation theory is valid. Perturbation 
theory holds for small go, or alternatively for 'large' p s . 
Consequently, direct perturbative results can be obtained 
in the renormalized-classical region with T -C p s - 

We discuss some important features of the perturba- 
tion theory here, with further details appearing in Ap- 
pendix^ For dimensions 1 < d < 2, there is long-range 
magnetic order for go < g c at T = 0, but rotation symme- 
try is restored at any T > 0. This singular phenomenon 
accounted for by a two-step integration procedure which 
has been discussed in detail in Sections 6.3.2 and 7.1.2 of 
Ref. 0: first we integrate out the modes with Matsub- 
ara frequency w„ ^ 0, and then subsequently perform a 
rotational average over the static modes by an exact pro- 
cedure. The first step is easily performed by a perturba- 
tion theory in which we assume that the local magnetic 
order is polarized along, say, the (0, 0, 1) direction. We 
obtain an expansion for the free energy in the presence of 
an applied magnetic field H a , which we assume has the 
value 

H a = (H ± ,0,H {[ ). (3.1) 
This expansion is discussed in some detail in Ap- 
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pendix IA II and yields the following expression for the 
free energy 



T = -T\vlZ 
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here J-q is the free energy in zero field. In (|3.2|) m has the 
apparent interpretation of the local magnetic moment of 
the impurity, while X±,\\ appear to be the transverse and 
longitudinal susceptibilities. However, it must be kept 
in mind that we are working in a T > regime where 
the magnetic order is ultimately averaged over and so to, 
X||,_l ar e merely intermediate quantities which arise in 
our computation, and do not have independent physical 
meaning. For go < g c the moment to is quantized exactly 
at the value to = S at T = 0, but corrections do appear 
at T > 0, as shown in Appendix IA II Following the 
method discussed in Section 6.3.2 of Ref. fl6l to the order 
in perturbation theory being considered here, the second 
step of rotational averaging over the directions of the 
local magnetization leads to the following expression for 
the physical magnetic susceptibility 
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Only the final quantity Ximp is a physical observable at 
T > 0. 

We can divide the contributions to the quantities in 
(|3.3|) to those arising from the bulk antifcrromagnet 
(which are proportional to its volume) and to those as- 
sociated with the impurity. First, for completeness, we 
recall results for the bulk susceptibilities, which are im- 
plicitly expressed per unit volume; there is no bulk con- 
tribution to the magnetic moment to. The results of bare 
perturbation theory for the bulk susceptibilities quanti- 
ties are given in (|A6I) and (|A7() . We re-express the results 
by replacing go by the physical p s ; these two quantities 
are related by>i& 
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In this manner, we obtain 
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Notice that both expressions have an ultraviolet diver- 
gence for d > 2, and so depend on the upper cutoff of the 
momentum integration. However, this divergence disap- 



pears in the physical bulk susceptibility 
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The last expression has been evaluated in d = 2, and 
we have re-inserted factors of c; this result has appeared 
earlier in the literatureiZ^. 

It is also interesting to see how ()3.6|l can also be ob- 
tained by the dimensionally regularized expressions in 
(IA6|) and (|A7|) . In dimensional regularization, the re- 
lationship Ij3.4[l becomes simply p s = X/go] substituting 
this into the integrals already evaluated in (|A6I) and (|A7|) 
we obtain 
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here T(s) is the Gamma function, and £(s) is the Rie- 
mann zeta function. Eqn. I|3.7|l agrees with Q3.6|l after 
using C(0) = -1/2. 

After subtracting out the bulk contributions to (|3.2(l . 
we are left with the impurity magnetization and suscep- 
tibilities. These can be computed by the same method 
as for the bulk susceptibilities. We will discuss the im- 
purity response to a uniform magnetic field at T > in 
Section IIII Al Section IIII Bl will consider the case of a 
local magnetic field applied only in the vicinity of the 
impurity site, while Section fill CI generalizes our results 
to a uniform magnetic field at T — 0. 



A. Impurity susceptibility at T > 

Evaluating first the frequency summations in 1|A9(I . 
then inserting (|3.4|l . we obtain the following result for 
the impurity magnetic moment, to, as an expansion in 
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Notice that while intermediate terms (like those in Q3.4[l l 
had a linear ultraviolet divergence in d = 2, the final 
expressions only have a logarithmic dependence upon 
the ultraviolet cutoff in d = 2. Unlike the case for the 
bulk susceptibility above, this divergence will not cancel 
against any other term. We can evaluate the integrals 
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in (|3.8[1 with a cutoff A in d 
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As we noticed above for the bulk susceptibility, the 
result (|3.9(1 can also be obtained in a somewhat simpler 
manner by the dimensionally regularized expressions in 
Appendix IA II Using the integrals already evaluated in 
[9l) we obtain 
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where 
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Finally, we take the limit d — ► 2 in (|3.11|) . This is found 
to be singular, as 771^2 both develop poles in (2 — d). In 
particular mi — > l/(27r(2 — d)) and m 2 — > — l/(47r 2 (2 — 
d)). As is conventional, we may identify the poles in 
(2 — d) with the logarithmic dependence upon the cutoff, 
and in this manner our earlier result 1(3. 9|) is seen to be 
perfectly consistent with l|3.10|l and (|3.11|) . 

We may proceed in a similar manner to an evaluation 
of the expressions for X-L,imp and XlLimp i n <|A10fl and 
(|A11|) . Here rather than using a momentum cutoff, we 
use the insights gained above to proceed with the simpler 
dimensional regularization method. Inserting the result- 
ing expressions into (|3.3|) we obtain the final result 
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As below (|3.11JI . upon taking the limit d — > 2, the expres- 
sions in (|3.13|) are seen to have simple poles in (2 — d) 
with xi - S ,2 /(3tt(2 - d)) and X 2 = -S ,2 /(6vr 2 (2 - d)); 
we replace the poles by ln(A/T) and thence obtain the 
result <|1.4|) announced in the introduction. We have also 
checked 1(1.41) directly in d = 2, by estimating the value of 
momentum integrals with a finite cutoff, using integrands 
similar to (|3.8|l obtained from Appendix IA II 

Close to the critical point, the expansion (|3.12(l implies 
that for general 1 < d < 2, the small y expansion of the 
scaling function <&(?/) in 1)1. 2f> has the structure 
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with the universal numbers xi,2 specified in l|3.13|l . The 
d — > 2 limit of xi,2 then leads directly to p. 3(1 . We are 
unable to obtain the values of the universal constants C 3 4 
because accurate results for the critical point are only 
possible for d close to 1, but in that case the logarithms 
of (|1.3|l are absent. 

A significant feature of these results for Xi mp 
is that while there is a (T/p s )\n(l/T) term, the 
(T/p s ) 2 ln 2 (l/T) terms have cancelled against each 
other; alternatively stated, the double pole in (2 — d) 
that is apparently present in l|3.13|l (associated with 
[r(l — d/2)] 2 ) turns out to have vanishing residue because 
C(0) = —1/2. This is an indication that this particular 
log singularity does not exponentiate upon inclusion of 
higher order terms, and is rather a consequence of Gold- 
stone spin-wave fluctuations, as opposed to a critical sin- 
gularity. 



B. Local susceptibility at T > 

We now consider the local susceptiblity, xioc which is 
the response to a field applied at the impurity site only. 
This is to be distinguished from the impurity suscepti- 
bility which is the response to a uniform field, after sub- 
tracting out the bulk contribution. The small go expan- 
sion for xioc is discussed in Appendix IA 21 The relation- 
ship (|3.3(l now generalizes to 
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and expressions for the terms on the r.h.s. appear in 
HA12fl . Evaluating the frequency summations and the 
momentum integrals with a cutoff A we obtain in d = 2 
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Here Sioc is a non-universal impurity moment which de- 
pends upon microscopic details like the local coupling 
constants and the precise location over which the field is 
applied. It is, in general, not equal to S, the moment 
which appears in Ximp- Similarly, C3.4 are non- universal 
numbers. Note however, that as in (|1.4|) . there is no term 
of order (T/p 2 ) ln 2 (l/T). 



C. Zero temperature response to an applied field 

The divergent impurity susceptibilities obtained above 
as T — > suggest that the response to a field will be 
singular at T = 0. 

At T — 0, the magnetic symmetry is broken for d > 1 
and small go, and so the quantities to, X||,impj and X-L,im P 
retain their separate physical identities and can be dis- 
tinguished experimentally. 

The calculation of the impurity response to a magnetic 
field at T = proceeds in a manner similar to that at 
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T > 0. The first crucial observation is that we now have 2 
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to all orders in g . This is a consequence of a 'gauge in- 
variance' of the action Z associated with the preserved 
symmetry of rotations about the z axis, and the transfor- 
mation I|A1|I . Explicitly, it is not difficult to check that 
upon converting the frequency summations to integrals in 
(|A9I) and I jAllj l. and evaluating the frequency integrals, 
the results in l|3.17[) hold to order g 2 , — this is a strong 
check on our computations. 

It remains to compute the transverse susceptibility 
X±.imp- Because of the broken spin rotation symmetry, 
this quantity is not protected by gauge invariance (the 
gauge symmetry is 'broken'), and it has non-zero contri- 
butions at each order in perturbation theory. However, 
certain terms in the perturbation theory have an infrared 
divergence for d < 2 in the presence of fully 0(3) symmet- 
ric Hamiltonian, and so we examine the full non-linear 
dependence of the impurity free energy on the applied 
field, J r imp(H±). The perturbative computation of Fhnp 
is described in Appendix I A 41 and from (1 A 16(1 we obtain 
for 1 < d < 2 
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This result, and the structure of the perturbation theory 
in Appcndix IA 41 suggest the follow universal scaling form 
for the critical behavior of T lrap {H j_) near the critical 
point: 
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the results here and in I imply that this scaling form holds 
for all 1 < d < 3. The results of I implicitly assumed 
an analytic dependence of on small H\, so that 

$jr(y) ~ y for small y. However, our computations here 
show that this analyticity holds only for d > 2, and that 
there is a leading non-analytic dependence with $jr(?/ — > 
0) ~ for 1 < d < 2. Precisely in d = 2, there is a 
pole in fi defined in (|3.18l) . and as in our discussion for 
m, this implies a logarithmic singularity: 



S 2 

My -> 0) = —y ln(l/y) + C T y 



(3.20) 



where Cjr is an unknown universal number. In the or- 
dered state in d = 2, well away from the critical point, 
we have from (|A16|I . and as in i|1.4fl 



^Fimp(Hx) — •? 7 imp(0) 



S 2 Hj fC 3Pi 

— In 

4irp s 

.90 < 9c 



0, (3.21) 



where C3 is a non-universal number. The logarithmic 
singularities in l|3.20|) and 13.21fl can be cut-off by spin- 
anisotropies in the underlying Hamiltonian, as has been 
illustrated in Appendix I A 41 



IV. RENORMALIZATION GROUP THEORY OF 
CRITICAL PROPERTIES 

There is already a well-established theoryi^ for the 
bulk phase transition at go — g c . Here we will show 
how this theory can be extended to the impurity corre- 
lations. This will be done with a single additional im- 
purity wavefunction renormalization constant Z' — from 
the perspective of boundary critical phenomena, this is 
a boundary renormalization factor at the impurity site 
x = 0. As noted earlier, the Berry phase m tSi m p is an 
exactly marginal perturbation to the bulk critical point: 
it is protected by its topological nature, and hence there 
is no additional coupling constant renormalization asso- 
ciated with the impurity spin. We will use this critical 
theory to obtain results for the impurity and local sus- 
ceptibilities at T > 0, and for the field dependence of the 
free energy at T = 0. 

First, let us recall the bulk renormalization theory from 
Ref. UM There is a field renormalization factor, Z, de- 
fined by 



0, (4.1) 

where NR a is the renormalizcd field. In the present quan- 
tum impurity context, this renormalization will be ad- 
equate at all spatial points away from the impurity, as 
has been indicated above. Second, Ref. [la has a coupling 
constant renormalization 



9o 



l-d 



5<2+l 



(4.2) 



where g is the renormalized dimensionless coupling con- 
stant, fi is a cut-off momentum scale, and 



Sd 



27T d / 2 



(27r) d r(d/2) 



(4.3) 



is a phase-space factor. To two-loop order and in the 
minimal subtraction scheme, these bulk renormalization 
constants are given byi^ 



Z = 1 + — 



3ff 2 



Zi = i + i + fL(i + e /2) 



where 



e = d-1. 

These constants give the beta function 

9 3 



(4.4) 



(4.5) 



(4.6) 
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which has a fixed point at g = g* which describes the 
bulk quantum critical point, with 



= e-e 2 +0(e 3 ) 



(4.7) 



Let us now turn to the impurity correlations. These 
require only an additional boundary wavefunction renor- 
malization which we define by 



N a (x = 0,r) = V^N^Jx = 0,r). 



(4.8) 



We discuss the computation of Z' in Appendix IA 31 where 
we find 



Z 3e 



(4.9) 



This renormalization constant implies that impurity spin 
correlations behave as 



(N a {x = 0,T)N a (x = 0,0)) 



where 



rj = e + r/ + [1(g) 



1 

TV 



dhxZ' 



9o 



(4.10) 



dg 



g=g 

e + f]-2n 2 S 2 e 3 + C(e 4 ). 



(4.11) 



Here r/ is the nearly-vanishing anomalous dimension of 
the bulk critical point which was mentioned in Sec- 
tion ITT1 — it controls the decay of N a correlations suffi- 
ciently far away from the impurity: 



(N a (x,T)N a (x,0)) 



1 



r t+ri 



.90 = 9c 



(4.12) 

The results in Appendix IA II can also be easily used 
to obtain an e expansion for the universal constant C± in 
(|1.5|) which determines the anomalous Curie response of 
Ximp at the critical point. We begin by substituting the 
renormalized coupling g defined by l|4.2(l into the dimen- 
siona ll y regu larize d expr ession for Ximp 

defined by lEHll . 

(|A9|) . (jAlOf) . and HA11|) . We expand the resulting ex- 
pression to order g 2 , and then expand the coefficient of 
each such term in powers of e. Consistency of the theory 
demands that poles in e cancel at this point, and this 
is indeed the case. Finally we substitute the fixed-point 
value g = g* in (|4.7I) into this expression. We find that 
all dependence upon fj, disappears at this point, which is 
strong evidence for the universality expressed by H1.2JI : 
the final expression then yields 



Ci 



3 



1 + 2e 



7T~ 

12 



0(e 3 



(4.13) 



As is also the case with the bulk exponents in the d = 1+e 
expansion, we do not expect the estimate of (|4.13|) to 
be accurate at e = 1. As discussed in I, we expect on 
physical grounds that S 2 /3 < d < S(S + l)/3. 



Next we consider the local susceptibility. As in I, this 
diverges near the critical point as 



Xloc 



T 



T 



(4.14) 



with $i oc a universal scaling function; the small argu- 
ment behavior of $i oc should be compatible with Ij3.16|l . 
while its infinite argument limit is a constant. By anal- 
ysis similar to that outlined in the previous paragraph, 
the expressions in Appendix IA 21 can be verified to be 
consistent with H4.14[) and the value of rf in (|4.11l) . 

Finally, we turn to the response to an applied field 
at T = 0, discussed earlier in Section Till CI and also in 
Appendix IA 41 At the critical point, any applied H± will 
induce long-range magnetic order in the bulk 19 for d > 1 . 
The scaling form H3.19JI nevertheless holds as p s — > 0, and 
we therefore obtain 



J~i mp (H±) — J-"imp(0) — —CjrH± 



go = g c (4.15) 



where = <&f(oo) is a universal number. This univer- 
sal number can be obtained directly from l|A16|) by the 
methods discussed above for Xi mp , and we obtain 



C 



e + 0(e 2 ) 



(4.16) 



V. CONCLUSIONS 



This paper has introduced the field theory (|2.3I) as a 
description of the low temperature properties of arbitrary 
static impurities in quantum antiferromagnets. The bulk 
fluctuations of the antiferromagnet are described by the 
familiar 0(3) quantum non- linear sigma model. Remark- 
ably, this venerable and strongly interacting field theory 
permits an exactly marginal perturbation, albeit on a 
'boundary', which has not been noticed before: this is 
the topological Berry phase of a spin S impurity. We 
have computed here the physical consequences of this 
marginal perturbation and so obtained a new description 
of the spin dynamics of the impurity. 

A preliminary comparison has been madeSi between 
our theoretical result (|1.4|) and the numerical results of 
Ref. 0. Reasonable agreement is found for the impurity 
susceptibility of a vacancy, and a more detailed compar- 
ison will appear later. 

Near the quantum critical point associated with the 
loss of long-range antiferromagnetic order, our results 
were obtained in an expansion in (d — 1). Key scaling 
features of these results were found to be in good accord 
with those obtained earlier in a (3 — d) expansion in I. 
In particular, right at the critical point, we confirmed 
the existence of a Curie 1/T impurity spin susceptibility 
but with an anomalous Curie constant not given by an 
integer or half-odd-integer spin. Our new result for the 
Curie constant is in (|4.13|) . Although the numerical es- 
timates for critical properties obtained from the (d — X) 
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expansion seem rather unreliable, this expansion never- 
theless provides convincing evidence for the existence of 
a strongly-coupled impurity fixed point with the physical 
properties discussed herein and in I. 
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APPENDIX A: DIAGRAMMATIC 
PERTURBATION THEORY 

We will consider perturbation theory in the presence of 
an applied uniform magnetic field H a under which (|2.3() 
is modified by 

d T N a — > d T N a — ie a p 1 HpN 1 

r-l/T 



co,k 



9o 



m p " ^ i m p S H a 



p-^x Q t dTN a {x = 0,r) (Al) 





This construction ensures that H a couples to a conserved 
total spin of the Hamiltonian. 

As in Refs. ll5il6l the perturbation theory in go is gen- 
erated by assuming that N a is locally polarized along a 
particular direction (say (0,0,1)), and by expanding in de- 
viations of N a about this direction. We do this here with 
the following parametrization in terms of a complex field 
ip, adapted from the Holstein-Primakoff representation: 

(A2) 

The advantage of the representation (|A2|) is that with 
the gauge choice 



A a (n) 



1 



l + n z 



-n y ,n x ,0) . 



(A3) 



the Berry phase takes the following simple exact form 

-5 (A4 > 

where the right-hand-side is to be evaluated at x = 0. 
Furthermore, the measure term in the functional integral 
also has the simple form 

J VNJ (Nl - 1) = 1 J (A5) 



co,k 





► 


• 


— >~-^ 


co,k 


► 




co,k 
— y 2 — 






■ 










ro 2 ,k 2 




_ k-| +K2~i~Kg 



»3' k 3 



2^2 



0)2+k2 

-S(-ico+H||) 
-2icoH||/g 

ifl>iH||/g 




co,k 



-CO1+CO2+CO3, 

> — □ — 



see 

co 3 ,k 3 caption 



co,k „ „ 
-> (H|| 2 -HL 2 /2)/g 



co,k 





& 2-i/2SH ± 5(co) 



FIG. 1: Propagator and vertices appearing in the calculation 
to the order needed. The weight of the fifth term is (&2 • &3 + 

L02L03 + kl ■ (k'2 + k 3 — kl) + U>i(lO'2 + LU3 - ui))/(4go)- 



(a) 



(b) 




FIG. 2: Diagrams for the bulk susceptibilities to order g$. 



The remaining terms in the action are obtained by insert- 
ing i|AHA2|) into l|2.3|l and expanding the results in pow- 
ers of tp — this yields a number of non-linearities which 
are analogous to those that appear in Ref. ITHl and these 
can be used to generate a Feynman graph expansion in a 
similar manner. We summarize the propagator and the 
vertices, to the order needed in our computation here, in 

pig. in 

First, we recall the results for the bulk response, in 
the absence of the impurity. The free energy is expanded 
as in (|3.2|) , and this leads to the diagrams in Fig [21 to 
order g®. There is no bulk linear dependence on H a to 
all orders in go-, and hence to = in the absence of the 



impurity. To quadratic order in H a we have the bulk 
susceptibilities (per unit volume) 

Xx,b = — + Hb) 
9q 



J2t) 



- T E 



2T d - 1 l 7 1 (27r)"- 2 C(2-d) 



(27r) d w 2 + k 2 



d-2. 



(A6) 



and 



X\\,b 

Gb) 



Gb) + Gb) 



2T V 

^ / (27r)<^ 2 4 



-1(2-1 



fc 2 



4T a - L J 1 (2TT) a - z ((2-d) 
d d k 



_4T V 

^ / (2^ 



_ 8T «-^ 2 ( 27 r) d - 2 C(2-(i) 



where 



- = r(o- 



d/2) 



(A7) 



(A8) 



(47r) d / 2 r(a) ■ 

As discussed in Section llTTl all intermediate Matsubara 
frequencies in all diagrams in this appendix are summed 
only over non-zero values; the integration over the zero 
Matsubara frequency modes leads to (|3.3|) . There are no 
infrared divergences in any graph (because of the summa- 
tion over non-zero Matsubara frequencies), while ultra- 
violet divergences appear in individual graphs for d > 1. 
We also list the expressions for the individual graphs 
obtained in the dimensional regularization method, ob- 
tained by analytic continuation from the d < 1 region — 
these will be useful in our renormalization group analysis. 
The dimcnsionally-regularized expressions ware obtained 
by first performing the momentum integrations, and the 
frequency summations are then naturally expressed in 
terms of the Riemann zeta function C( s )- There are 
also many sensitive cancellations in the ultraviolet diver- 
gences of the various graphs considered in this appendix, 
and these will appear as cancellation of poles in the di- 
mensionally regularized expressions. In Section II I II we 
have also considered the expressions of this appendix di- 
rectly in d = 2 without dimensional regularization, and 
these results illustrate the cancellation of ultraviolet di- 
vergences upon expression of the results in terms of phys- 
ical observables. 

The application of the perturbation theory towards 
computation of physical properties of the impurity in dif- 
ferent regimes will be presented in separate subsections 
below. 



1. Impurity susceptibility at T > 




FIG. 3: Diagrams contributing to m to order g^. 



The diagrams for the perturbative expressions for the 
impurity contributions to the quantities in (|3.2I) are 
shown in Figs 13151 Now there is a contribution to linear 
order in H\\, and Fig yields the following expressions 
for m: 



Eh) 



Gb) 



Qb) 



s + Hi) + ©) + ©) + ($) 

~2Sg T d - l J 1 {2n) d - 2 C{2 - d) 

45. go T^ 1 1 7 2 (27r) d - 2 C(2- ( i) 

d d h w 2 



-2Sg 2 T 2 



E 



(27T)" {Ul+klf 



E 



d d k 2 



(2n) d el + kl 



m = ^g 2 T 



&SglT 2d - 2 J x J 2 [(2^) d - 2 C(2 - d)] 
d d k ujI 



E 



(2ir) d [u>l + k 2 ) 2 



l6Sg 2 T 2d - 2 J 2 [(27r) d - 2 C(2 - d)] 2 . (A9) 



We first address the computation of the impurity mag- 
netic susceptibility, Ximp a t non-zero temperatures. 



Similarly, for x±,imp we only have the diagram in Fig0| 
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FIG. 4: Diagram contributing to x±.im P to order g$ at T > 0. 




(a) 



- o 

- o 



(e) 



FIG. 5: Diagrams contributing to %\\ imp to order g$. 




which we obtain: 



X|| ,imp 

Eh) 



Eh) 



Eh) 



Eh) + Eh) + Eh) + Eh) + &) 

1 

x R+l 

2S* 2 .ggT 2d - 3 J 1 2 (27r) 2d - 4 C(4 - 2d) 

45 2 . 9o 2 T x: r ddki ddk2 



(2n) d (2n) d {uj 2 n + fc 2 ) 2 
1 



K + fc 2 ) 2 

8S 2 ^T 2d - 3 1 7 2 2 (27r) 2d - 4 C(4 - 2d) 



8S 2 g 2 T £ 



d d fci d d fc 2 



^7, 



(2n) d (2ir) d (wg + fc 2 ) 3 
1 



(^ + fc 2 2 ) 



16^^T^- J J 3 J' 1 (27r)^- 4 C(4 - 2d) 
d d /fci d rf fc 2 cj 2 



8 o 2 a 2 T V /• 



X 



) d (^ 2 + fc 2 ) 2 



(^ 2 +fc 2 2 ) 



\QS z g z Q T M ~ A JiJ 2 (27r)^" 4 C(4 - 2d) 
d d fci d d fc 2 w 2 



2^ g Q ± ^ J {27r )d {2 n 



X 



{2n) d (2ir) d (w* + k 2 ) 2 
1 



K + fc 2 2 ) 

~AS 2 g 2 T 2d - 3 J 1 J 2 (2ir) 2d - 4 ((4 - 2d) (All) 



which yields: 



X-L,i m P 







d d fc! d d fc 2 



> d (2^) d (co 2 + /fc 2 ) 2 



(u> 2 +fc 2 ) 

25 ,2 ff 2 T 2<i - 3 l 7i.7 2 (27r) 2 ' 1 - 4 C(4 - 2d).(A10) 



2. Local susceptibility at T > 

The response to a field applied only near the impu- 
rity site can be computed as in Appendix IA 21 Only the 
graphs in Fig [3h and Eh now contribute, and we have 
therefore 

rn-ioc = S + Eh) 
xii.ioc = Eh) 



X-L,loe = 



(A12) 



where the values of the respective graphs are as specified 
in and l|A"TT]) . 



3. Spin correlations at T — 



Finally, for XlLimpj we have the diagrams in FigEl from 



The methods above can also be extended to obtain im- 
purity spin correlations at T = and g = g c . As long as 
we restrict ourselves to rotationally invariant correlation 
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FIG. 6: Lowest order diagram contributing to the correlator 
in 1 A 1311 which depends upon the presence of the impurity. 
The X'a denote external sources for the fields. 



functions, direct perturbation theory in go is free of in- 
frared divergences. For the impurity spin correlation in 
(|4.1U|) , the first corrections which depend upon the pres- 
ence of the impurity do not appear until order g^: these 
arise from the graphs shown in Fig H3 and lead to the 
following expression: 

(N a (x = 0,r)N a (x = 0,0)) = 1 + 
duj d d h d d k 2 d d k 3 



2tt (2n) d (27r) d (27r) d 
uj 2 (1 — cos(wr)) 
{uj 2 + k 2 )(uj 2 + k 2 )(uj 2 + 



(A13) 



Here the ellipses denote numerous lower-order terms 
which do not depend upon the presence of the impurity 
and hence are the same at x — and x ^ 0; the first 
term which breaks translational invariance is shown in 
(|A13|) . The integrals in (|A13|I can be easily evaluated in 
dimensional regularization, and the second term in (1 A 13(1 
equals 

2glS 2 [T{{2 - d)/2)fY{M - 3) cos(37r(d - l)/2) 
(47r) 3rf / 2 7rr 3d - 3 

(A14) 

Picking out the pole in e in (|A14|) , we immediately obtain 
(Oil . 



4. Response to a field at T = 

As discussed in Section IIII CI we need the impurity 
contribution to the free energy in the presence of an ap- 
plied transverse magnetic field H±, J-i mp (H±). We will 
see that the response is singular as H± — ► 0. The sin- 
gularity can be cutoff by an easy-axis spin anisotropy, 
and for completeness, we perform the computation in the 
presence of such an anisotropy. So we modify the action 
by 

S b [A^] S b [Ay -j J d d x J J drN 2 . (A15) 

Because we are now computing the free energy to all 
orders in the applied field, the Feynman graph expansion 
is quite tedious, and we will be satisfied by obtaining the 
result only to order g . The computation is done most 
simply using the Cartesian components ~ (tp+tp*, —i{ip— 
ip*)), and to leading order in go, only the graph shown in 



FIG. 7: Diagram contributing to the free energy J r (H±) to 
order go- at T = 0. The propagator represents the x compo- 
nent of N a and includes the open square vertex in FigQto all 
orders in H±, along with the easy-axis anisotropy in 1A15L 
it equals go/(u 2 + k 2 + D + H\). 



FigEJcontributes. Note that this diagram did not appear 
in the computation at T > because of the restriction 
there to summation over non-zero Matsubara frequencies 
and the delta function in frequency associated with the 
vertex in Fig [7] the leading term in x^imp was of order 
t/o at T > 0. From the diagram in Fig [7] we obtain 



-?"imp (H±) — — 



d d k 



goHjS 2 

2 J {2Tr) d k 2 

g H 2 ± s 2 



1 



D + H\ 
J^D + Hl)^- 2 ^ 2 . (A16) 



This graph has a log singularity in d = 2. The same 
logarithm appeared in a different manner in the T > 
computation: it was present in Fig [3^. Ultimately it is 
only Ximp that is physically measurable at T > 0, and it 
is clear now that the logarithm appears in different places 
depending upon the different organizations of perturba- 
tion theory at T = and T > 0. 



APPENDIX B: PERTURBATION THEORY FOR 
GENERAL 7 

The computations elsewhere in this paper have been 
limited to the case in which the coupling between the 
impurity spin and the bulk antiferromagnetic spin fluc- 
tuations, 7, has effectively been sent to infinity. We have 
argued that this limit appears naturally in the vicinity of 
the quantum critical point. This appendix will consider 
the general 7 case, and consider the extent to which the 
low T properties away from the critical point are inde- 
pendent of the value of 7. 

We shall be concerned here with the partition function 



3L 



VN a (x,r)S (Nl - 1) Vn a 6(n 2 a - 1) 
x exp (-Sb [AJ - Sim P , 7 ) 



l/T 



dr 



iSA a (n) 



dn a (r) 
dr 



jSN a (x = 0,r)n a (r) 



(Bl) 



where 5b [N a ] is as in (|2.3|) . and A a (n) is defined by (|2.2|l . 
In principle, it is possible to generate an expansion in 
powers of go, with each term containing its exact de- 
pendence on 7 and S; this requires an exact treatment 
of the impurity spin fluctuations, and this can be done 
by the method described in Appendix C of I. Here, we 
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shall use the method described above in Section 1X1 with 
a parametrization similar to (|A2|I applied also to n a . By 
this method, it is not difficult to obtain results order- 
by-order in go, dropping only diagrams with a 'tadpole' 
factor of the impurity spin propagator (i.e. with a sim- 
ple closed loop of the impurity spin propagator) — these 
are easily seen to have a prefactor of e~ 7 / T (we assume, 
without loss of generality, that 7 > 0). 

We now present results to order g for the impurity 
spin susceptibility at T > 0, computed above in Ap- 
pendix We will omit all details and merely present 
final results to leading order in go. Dropping terms with 
a pre-factor of e _7 / T , we found 



1 



2lLU n 



[ojI + k 2 )(-iuj n + 7 ) 2 (a* + k 2 ) 2 (-iuj n + 7) 

d d k 



A 1 1, imp 



= 2 1 SgoTY J J 



1 



(27r) d 



7 



K + fc2)3 ( W 2+fe2)2 ( W 2 +k 2^_ i0Jn + 1 y 
, 7 . 2i l^n 

{ujI + k 2 ) 2 (-iuj n + 7) (w 2 + k 2 ) 2 (-iuj n + 7 ) 2 

47^n 

(ojI + k 2 y{-iuj n + 1 ) 
X_L,im P = jSg T J 



7 



d d k 

(2n) d 



1 



(col + k 2 ) 2 



(uj 2 l + k 2 ) 2 (~iuj n + 1 ) 



(B2) 



It is now easy to check that the 7 — > 00 limit of these 
expressions is finite, and indeed agrees precisely with the 
order go results for Ximp obtained in Section IIII Al and 
Appendix IA II this is a non-trivial check of our compu- 
tations. Evaluation of the frequency summations in 1B2I) 
is a tedious but straightforward exercise. After this, we 
combine the results using (|3.3[l , and evaluate the momen- 
tum integrals at low T as in Section IlIII while keeping 7 
finite; in the limit of T <C 7, A we obtain in d = 2 



Ximp — 



T 

TTPs 



(B3) 



Notice that the co-efficient of the (l/p s ) ln(l/T) is in- 
dependent of 7, and that it agrees with Q1.4|) . Also, at 
finite 7, the Tln(l/T) term does acquire a non- universal 
7-dependent correction. 



APPENDIX C: LOW TEMPERATURE 
PROPERTIES IN d = 3 

This appendix briefly describes the extension of our 
results to d = 3. The bulk quantum critical point in 
d = 3 does not satisfy strong scaling properties, and so 
we will not consider it here. We will focus only on the 
low T properties within the magnetically ordered state, 
well away from any quantum critical point. 

Magnetic long-range order is present for a finite range 
of T > 0, and so the magnetic response remains 
anisotropic as T — > 0. The quantities m, X-L.imp, XlLimp 
retain their separate physical identities, and can be mea- 
sured separately. 

The low T expansions for m, X-L,imp> XlLimp are ob- 
tained as in Appendix ^ Indeed, now we need not sep- 
arate the oj n = and the u> n ^ modes as there is long- 
range order for T > 0: the expressions in Appendix lAl 
can therefore be used here, after converting all frequency 
summations to run over both zero and non-zero values of 
Matsubara frequencies. In this manner, (|3.8(l is modified 
to 



S 



1 



T 

Ps 



d 3 k 



1 



(2tt) 3 AT 2 sinh 2 (fc/(2T)) 
T f d 3 p ( 1 
7s J \pT(eP/ T - 1) 

1 



AT 2 sinh 2 (p/(2T)) 



(CI) 



Evaluating the momentum integrations, and re- inserting 
factors of c, we obtain 



S 1 



T 2 



T 



6cp s 72c 2 p 2 J ^ ^ 

Interestingly, the expression l|C2(l can also be obtained 
simply by setting d = 3 in (|3.10() . For the finite 7 case, 
discussed in Appendix ]B\ the T 2 term above remains 
unchanged, while the T 4 term does acquire 7-dependent 
corrections. 

The results for Xy^mp and Xx.imp now follow from 
(|A"TTJl and (JMUJ). Setting d = 3 in the dimensionally 
regularized expressions here, we find that the co-efficient 
of the universal T 3 term vanishes for both quantities. 
However, there are non-universal T 2 corrections for both 
XlLimp and X-L,imp; and these have to be estimated di- 
rectly from the expressions in (| All|> and IjAlOfl : the fre- 
quency summations have to be evaluated first (including 
the zero Matsubara frequencies), and then the momen- 
tum integrations have to be evaluated with a finite cutoff. 

Similar techniques apply to the response to a local field 
discussed in Appendix IA 21 Now we obtain the universal 
correction 



m\oc{T) - mioc(0) = - 



ST 2 
\2cp s 



(C3) 



along with non- universal T 2 corrections to XJ-.ioc and 
XII, loc- The factor of 2 difference between the T 2 terms 
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in (|C3|) and i|C2(l is an interesting characteristic of the 
theory. 



APPENDIX D: COMMENT ON SUSHKOV'S 
COMPUTATION 

It has been claimed by Sushko^ that the Curie con- 
stant remains C\ = 1/4 for a S = 1/2 impurity at the 
quantum critical point of an antiferromagnet. Here we 
show using the model of his paper that there is a pertur- 
bative correction to the impurity susceptibility, and that 
this implies an anomalous Curie constant. Of course, 
the possibility remains open that the (3 — d) and (d — 1) 
expansions both fail in d = 2 near the critical point, 
but reasons for such a possible failure do not appear in 
Sushkov's arguments. 

Sushkov models the bulk spin fluctuations at the quan- 
tum critical point using a S = 1 boson t a , as in Ref. l2ll 
These bosons are coupled to the external magnetic field 
H (assumed oriented along the z axis) and to the impu- 
rity moment S a - This gives us the model considered by 
Sushkov: 

Ti. = TIq + Ti-i 

Ho = (£(k)-mtf)4(k)i m (k)-if4 

k m=0,±l 

Hx = \4>J a (Dl) 

where k is the momentum of the t a bosons with energy 
e(k), and 

Q = ^-i=(i a (k)+tt( k )) ( D2 ) 
k V 2e ( k ) 

and 

t x = {t x +t-i)/y/2 

ty = - ^/V^ 

t z = t (D3) 

It can be checked that H couples to the total spin, which 
commutes with the Hamiltonian. 

Now we compute the free energy, J 7 , in a power series 
in A in arbitrary H. To second order in A, this is done 



by the familiar formula 

T = -TlnTre~ Wo/T 

- T L dT L dTi — — (D4) 

where 

Hi(t) = e HaT Hie- n ° T (D5) 

Everything in (|D4(I and (|D5I) can be evaluated analyti- 
cally by simple means, and then we can perform the in- 
tegrals over r and t\ - this was done using the computer 
program Mathematica for arbitrary H and S, without us- 
ing any diagrammatic perturbation theory. Finally, we 
can expand the result in powers of H and obtain for the 
impurity susceptibility 

_ S(S+l) 2X 2 S(S+1) ^ e £ ( k >/ T 
Ximp ~ 3T + 3T 2 ^ e(k) 2 (e e ( k )/ T - 1) 2 

(D6) 

This result agrees precisely with that obtained using a 
diagrammatic approach in I. It disagrees with that of 
Sushkov, who did not obtain any correction to the first 
free moment term — he does not appear to have consid- 
ered the cross-correlation between the bulk magnetiza- 
tion of the t a and the impurity magnetization. Note 
that this disagreement appears already at the level of 
bare perturbation theory, and does not involve any of the 
subtleties associated with approaching the scaling limit 
at the critical point in the e or 1/N expansions. 

In the quantum disordered regime above the paramag- 
netic phase, we can model e(k) = \J c 2 k 2 + Aj,, where 
A T ~ A > is the spin gapSiSi; here l|U6(l predicts a 
contribution of order e _A / T to the susceptibility, and so 
the moment is indeed precisely quantized at S. 

However, the quantum critical region^ we have At ~ 
T, and then (|D6(1 yields a contribution to the suscep- 
tibility of order X 2 T d ~ 4 . For d < 3, the dimensionless 
combination A 2 T d ~ 3 approaches a universal value at the 
fixed point, and a universal irrational correction to the 
Curie term applies, as shown in much detail in I. 
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